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Abstract 

A set of coupled time-dependent Ginzburg-Landau equations (TDGL) for 
superconductors of mixed d- and s-wave symmetry are derived microscopically 
from the Gor'kov equations by using the analytical continuation technique. 
The scattering effects due to impurities with both nonmagnetic and magnetic 
interactions are considered. We find that the d- and s-wave components of the 
order parameter can have very different relaxation times in the presence of 
nonmagnetic impurities. This result is contrary to a set of phenomenologically 
proposed TDGL equations and thus may lead to new physics in the dynamics 
of flux motion. 
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I. INTRODUCTION 



There are growing experimental evidences to suggest that high-T c superconductors have 
a dominant rf a .2„„2-wave pairing symmetry.0 Based on symmetry considerations, Voloviki 
argued that an s-wave component of the order parameter should be generated near the 
core region of a vortex in a <i-wave superconductor. This conclusion was later confirmed by 
a numerical calculation! and by studying a set of microscopically-derived two-component 
Ginzburg-Landau (GL) equations.^ 

In view of the enormous success of the GL theory for describing the equilibrium proper- 
ties of superconductors near T c , it is natural to generalize it to time- dependent situations. 
This generalization has become particularly desirable, since a set of phenomenological time- 
dependent Ginzburg-Landau (TDGL) equations for coupled s- and d-wave superconducting 
order parameters has been recently proposed,! and used to investigate the dynamics of vor- 
tices in high-T c superconductors. One would very much like to know how valid is such an 
approach. 

It is well known, however, that TDGL equations are not as universal in form as its 
time-independent variety, but can be dependent strongly on whether the system is gapful or 
gapless, and in the later case, whether a strong or weak gaplessness condition is assumed. 
The simplest set of TDGL equations for conventional s-wave superconductors was proposed 
phenomenologically by Schmid,! and subsequently derived microscopically by Gor'kov and 
EliashbergEI under the assumption of a strong gaplessness condition (i.e., t s T c <C 1, where t s 
is the spin-flip lifetime and T c is the transition temperature). This set of equations has been 
used in the past to study the vortex dynamics in conventional superconductors.! EliashbergP 
has later derived a more complex set of TDGL equations for low-T c s-wave superconductors 
assuming only the weak gaplessness condition [i.e., t s A <C 1, where A (T) is the equilibrium 
value of the (s-wave) order parameter in the absence of fields]. It has been used to study 
flux-flow resistivity^ and the transport entropy of vortices.EJ Even more complex sets of 
TDGL equations have been derived subsequently assuming only the dirtv limit condition 
(t]T c <C 1, where T\ is the total scattering lifetime) and (1 — T/T c ) <C 1,E3 so the system 
need no longer be gapless. However, this set of equations is so complex that it has not yet 
been extensively used. 

In this work, we shall derive microscopically a set of coupled TDGL equations for super- 
conductors with mixed d- and s-wave pairing symmetry based on the approach of Gor'kov 
and Eliashbergl'l in the presence of impurities with both spin-flip and non-spin-flip interac- 
tions and assuming only weak gaplessness conditions for both waves (i.e., TiA d0 <C 1, and 
r s A s0 "C 1.) The primary objective of this derivation is to establish a reasonably reliable 
set of equations governing the dynamics of coupled d- and s-wave order parameters which 
are hopefully valid for describing the dynamic properties of high-T c superconductors. 

The outline of this paper is as follows: In Sec. II, the TDGL equations for the order 
parameters are derived. The expressions for current and charge density are presented in Sec. 
III. Finally, discussions and summary are given in Sec. IV. 
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II. TIME-DEPENDENT GINZBURG-LANDAU EQUATIONS FOR THE ORDER 

PARAMETERS 



We begin with the Gor'kov equations!! 

[-£ - / l (xr)]G^(xr,xV) - £/ Q7 (x)G 7/3 (xr,x'r') 

+ / dx"A a7 (xr°+, x"r)F 7 yx"r, x'r') = <5(x - x')5(r - r')<W , (2.1a) 
[£ - /^(xr^^xV) - [/ 7Q (x)F 7 y X r,xV) 

- J rfx"At 7 (xr°+, x"t)G 7/J (x"t, x'r') = . (2.1b) 

Here repeated spin indices mean summing over these indices. In these equations, 

h(xr) = ^ 1 n - e<p xr - /i , 2.2 

2m 

is the single-electron (— e) Hamiltonian with A(xr), </?(xr), and /i denoting the vector, scalar, 
and chemical potentials. (We have assumed h = c = 1.) By assuming zero-range interactions 
bewteen electrons and impurities, the impurity scattering potential can be written as 

C/ a/3 (x) = Y^PiSap + U 2 (S t • ^)]5(x - Rj) , (2.3) 
is/ 2 

where I denotes the set of impurity sites, cr is made of the Pauli spin matrices, Sj is the 
spin carried by an impurity at Rj. U\ and U2 are the non-spin- flip and spin-flip interaction 
strengths, respectively. By definition, the order parameter in real coordinate and imaginary 
time space is 

A^xr, x'r) = V (x - x')Flp(xr 0+ , x'r) , (2.4) 

where — V(x — x') is the effective pairing interaction between electrons. Because of the 
spatial and temporal non- uniformity, the Green function G aj g(xr, xV) and F^x^xV) 
are not the functions of coordinate and time differences. When expressed in the imaginary 
frequency space after the Fourier transform, they depend on two frequency variables. For the 
spatial coordinate dependence, as treated in the static case,@ we express these two functions 
in terms of the center-of-mass coordinate R = (x + x')/2 and the relative momentum after 
a Fourier transform with respect to the relative coordinate r = x — x'. Thus Eq. ( [2.4]) can 
be rewritten as: 

r r/k' 

A^R, k; u) = T £ J j^V(k - k')F^(R, k'; e, e - u) , (2.5) 

where k is the relative momentum, uo = 2imrT and e = i(2n' + l)7rT with integers n and n', 
— V(k — k') is the pairing interaction in the momentum space, ^'(R, k; e, e') is the Fourier 
transform of F^(xr, x'r'). To relate to high-T c superconductors, we have assumed that 
the system under consideration is two dimensional. For the spin-singlet pairing, the order 
parameter is given in the spin space as A*^ = A* g a p, where 

9a0 = f _°! J ) • (2-6) 

V / a/3 
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To obtain the TDGL equations for superconductors of a mixed d- and s-wave symmetry, we 
make the following ansatz for the pairing interaction and the order parameter 



V(k - k') =V. + V d (kl - - K") , (2.7) 



A*(R, k; u) = A:(R; u) + A*(R; uj)(k 2 x - k 2 y ) , (2.8) 

where Vd and V s are positive so that both the d- and s-channel interactions are attrac- 
tive. The (i-channel attractive interaction could originate from the antiferromagnetic spin 
fluctuation, whereas the s-channel attractive interaction might arise from phonon mediation. 
Introducing the Green function G° of the normal metal, which satisfies the equation 



d_ 

"dr 



h(xr)]G%(xT, xV ) - f4 7 (x)Gyxr, xV) = 5(x - x')5(r - r')5 aP , (2.9) 



Eq. (|2.1|) may be converted to a set of coupled integral equations: 

G^(xr,xV) = G° /3 (xr,xV) 

-J dx 1 dx 2 dT 1 Gl^xT,x 1 T 1 )A^ u (x 1 Ti,X2Ti)Fl l3 (x 2 T 1 ,x'T') , (2.10) 

F^(xt,xV) = J rfxirfx 2 rfriG° a (xiri,xr)A* ;y (xiri,X2ri)G i , /3 (x2ri,x / r / ) . (2.11) 
Also note that the normal-state Green function can in turn be written as an integral equation 
Gyxr.xV) = GyxT.xV) + J dx"dr"Gl ( (xT,x"r") 

,,r e ^"( x T ) ' P x " I 'I "\l/"<0 I II II I i\ /<-> io\ 



with G° as the normal-state single-particle Green function in the absence of the electro- 
magnetic field but including the effect due to impurity scatterings. To write down the 
above integral equation, the squared term of the vector potential has been neglected and 
the Coulomb gauge is chosen. 

A. Analytical continuation 

To incorporate the time dependence of physical quantities, we use the analytical continu- 
ation technique discussed in Refs.BH to transform imaginary frequencies into real frequencies. 
The procedure is as follows: (i) In Eq. ( |2.5|) , each term of the summation over the imaginary 
frequency e can be regarded as the residue of an integral along the contour around the point 
z = e so that we have the transformation ~~ > 4^ §c dzt&rih. Associated with this 
transformation, all involved e are replaced with z. For example, 

Tj2G°(-e)G°(e-iu') = — i dzt&nh^-G°{-z)G°{z - u') , (2.13) 
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where the spatial and spin variables have been suppressed for simplicity, (ii) Deform the 
contour integral around z into the straight line integrals along z = e±iO + , z = e+u'±iO + , . . . , 
where e G (—00, 00) is the real integral variable and + is infinitesimal. So far, u', ... are 
still imaginary frequency 2nniT and we take n > since we will perform the analytical 
continuation from the upper half-plane. As a consequence, each Green function G°(z) with 
energy variable coincident with this line is decomposed into G°^- ) (e)— G°^(e), where G 0( - R ' A ^ 
are the retarded and advanced Green functions, respectively. The minus sign before G ^ 
comes from changing the direction of integration. The other Green functions are mapped to 
the retarded or advanced Green function depending on their energy variable. Then Eq. Q2.13 ) 
becomes 

£ JZo detanh #{[G°^(- e ) - G^ A \-e)\G^ A ) {e - uj') 

+ G°( R \-(e + uj'))[G°^(e) - G°^(e)}} . (2.14) 

(iii) Since all uj', . . . , lie in the upper half-plane, the Fourier transform of the expan- 
sions ( ^.lUp and (2.11) are analytical in each of these variables. Therefore, we can implement 



the analytical continuation by simply replacing all uj 1 with uj % + i0 . Simultaneously the dis- 
crete summation TY^ U ' is replaced by a continuous integral (27r) _1 / did and the Kronecker 
delta T~ x b w ^i by the Dirac delta {2ti)8{uj — uj'). Finally, we get 

h S-oc de \~ tanh ^fG°^(-e)G°^(e - d) + tanh ^fG^ R \-e)G^ R \e - uj')} 

+ 4^/"oo^G°^(-e)[tanh^-tanh^]G°( A )( e -^) . (2.15) 

The first two terms consist of the only advanced and the only retarded Green functions 
(Following Ref.B, we shall refer to them as the normal part). The rest part has those terms 
involving the product of retarded and advanced Green functions, in which a change from a 
retarded to an advanced Green function occurs in only one place (We shall refer to them as 
the anomalous parts). After obtaining the results in terms of real frequency, we can perform 
the inverse Fourier transform to represent them in a real time. 
Using these rules, we can obtain the following expression 

T £ Ft (e, e - d) - ±- J de tanh ± [F*W (e + u,, e) - (e, e - u,)] 
+F^ R \e,e 1 )A^ A \e 1 -d,e-u)] 

+ [G + ™(e, e 1 )F^ A \e 1 -d,e-d) + F^ R \e, e l )G^ A \e 1 -d,e- uj)] 
x [eA^ ■ p/m - ecpuf] } [tanh ^ - tanh ^ ^ ] , (2.16) 

where, G~^ R ' A ^ and F^ R ' A ^ are formally defined by 

G"(p, p - k; e, e - uj) = G°-(p; e) + G°-(p; e)A £J ,(k / )G°-(p - k'; e - w')A^(k") 

xG°-(p - k' - k"; e-uj' - d') + . . . , (2.17) 

and 
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Ft(p, p - k; e, e - u) = -{G°+(p; e)A:,(k')G°-(p - k'; e - a/) + . . .} , 



(2.18) 



in which the substitution e — > e ± i0 + for the retarded (advanced) Green function should be 
made. Here G 0± (p; e) = [e±^ p ] _1 and all uj 1 are real. The functions G + and F are obtained 
from G~ and F 1 " by changing the sign of £ in G 0± . Note that G + and F are introduced only 
for simplicity of notation. 



B. Normal part 

The normal part of A* can be written as 
dk' 

( 27F )" e„>0 e„<0 

(2.19) 



A^(R, k; W ) = / 7^^(k " k')[T £ F™(R, k'; e n + U , e n ) + T £ ^(R, k " e n - c 



The evaluation of the normal part can be done by expanding the expressions in powers 
of the order parameter. We write for F^ R,A) and G^ R,A) up to terms of the third and second 
order in A, respectively, so that 

Flf^ (x, x'; e, e - u>) = (x, x'; e, e - u) + F^f (x, x'; e, e - u>) 

+^& A) (x,x';6,6-^), (2.20) 

where 

F%$\x,x!;e,e-u>) = J d^ 2 G°^ A \^, x; -e)e'- A ^<^ A*^, x 2 ; u") 

xGffi A) (x 2 ,x!-, e - u ) e -*K>»(x*)<**-*') ; (2.21) 

F^f (x,x';e,e-c) = | dXidxac^GjL^Cxx.x; -e)A^( Xl , x 2 ; u/)G#^(x 2 , x 3 ; e - a/) 

x [-e</v (x 3 )]G°ff ' A) (x 3 , x'; e - uj) 

+ J rfx 1 rfx 2 rfx 3 G^(x 1 ,x 3 ;-(e-cc;0)[-e^(x3)]G^Hx3,x; -e) 
xA^x.x^'^^^x.x'^-w) , (2.22) 

F ii R af(^ x '; e> e - cj) = - y dx 1 dx 2 dx 3 tix4tZx5dx 6 G ! 2f ,A) (xi, x; -e)A* i/ (xi,x 2 ; u/) 

xG°?' A) (x 2 , x 3 ; 6 - c/)A A(7 (x 3 , x 4 ; cu")G^(x 5 , x 4 ; -(c - a/ - a/')) 

x A; r (x 5 , x 6 ; OGjr^' x '' e " w ) • ( 2 - 23 ) 

Here the summation over the imaginary frequency w 1 with the constraint J2i ^ = w is i m_ 
plied. To write down the above expression for F^ R ' A ^ , we have expanded G° to the first order 
in the scalar potential — eip, and separate this expansion term out explicitly. As far as the 
dependence of G° on the magnetic field is concerned, the quasiclassical phase approximation 
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can be used to write it in the form G°(x, x'; e) exp[— ieA^ • (x — x')]. Accordingly, the gap 
function can also be written as a sum of three parts 

A^(R, k; u) = A*f Q/3 (R, k, u) + A%%(R, k; u) + AJ>%(R, k; u>) . (2.24) 

The remaining task involves the evaluation of the average over an ensemble of randomly 
distributed configurations. As an approximation, A* is regarded as very nearly independent 
of impurity configurations. We assume that impurities density n« are randomly distributed 
and their spins are arbitrarily oriented so that there is no correlation among them.0 Using 
the Born approximation, we can show that the impurity averaged zero-field normal-state 
Green function takes the following form 

(Gyx, x'; e n )) = / dke^-^ , (2.25) 



where (. . .) denotes the average over the impurity configuration, £k = ^- — ft is the kinetic 
energy, and r\\ — 1 + (2r 1 |e n |) _1 with the scattering time T\ given by 

- = 2v™ J iV(0)[|f/ 1 | 2 + -S{S + \)\U 2 \ 2 \ . (2.26) 

Ti 4 

Here rti is the impurity density and N(0) is the density of states at the Fermi surface per spin. 
The evaluation of the product of Green functions can be conveniently performed based on 
the diagrammatic ruleBEl If there are two Green functions connected by an s-wave order 
parameter or s-channel two-body interaction, these two Green functions might be called 
directly connected and we should attach a vertex renormalization factor: 

^)(e) = [l- fl (^)(0)]- 1 

g -J- i 

« 7±X T ( 2e ±^) 2 ' (2 ' 27) 

where 

a<R, (°» = 2^/ wf™^ - • (2.28) 

and 

= « / Wf- G ^ - £)G ° W,(P - £ + iw) ' (2 - 29) 

Here 

- = 27rn J iV(0)[|t/ 1 | 2 - ]s(S + 1)\U 2 \ 2 ] , (2.30) 
r 2 4 

and the spin-flip scattering rate is defined as 2T" 1 = rf — r 2 _1 . If the two Green functions 
are connected by a (i-wave order parameter or d-channel two-body interaction, they might 
be called not directly connected and we have no vertex correction. For the average of the 
product of more than two Green functions, an impurity line can also appear across the box 
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of a diagram. Because of their zero contribution, diagrams with more than one impurity 
line across the box should not be included. This impurity averaging technique was first used 
by Abrikosov and GorkovB for conventional s-wave superconductors. Recent experimental 
measurements by Bernhard et a/.0 on various types of YBa2(Cui_ x Zn ;c )307„5 samples have 
shown that the depression of T c by Zn doping can be fitted well with the Abrikosov- Gorkov 
theory applied to the (i-wave superconductivity. 

Now we give a derivation for the gap function from TJ2e n >o F^ R \e n + to, e n ). The con- 
tribution from TJ2e n <o F^ R K e n, € n — id) can be obtained by merely changing all explicit 
i to —i and id to — u, which gives the same result. In addition, one can easily see that 
Aj^JR, k; u) = since the contribution from the two terms given by Eq. ( |2.22| ) cancelled 
with each other. Therefore, we obtain 



A^(R,k;c) =T £ / T^(k-k') / dxe*"J J dR'dr' 



r' r _ „nrm ,_, r'_ r 



x (G°£\B! + - , R + -; -(e n + id))g x ,G% R \R> e»)> 
xe^'- R )- n /^e lk "- r 'A:(R,k"), (2.31) 



and 



A*^Jf (R,k;w) =T J 7^¥ V ( k - k ') / dre ~ tk ' r J dR'dr'dR 1 dr 1 dR 2 dr 2 



Cn>0' 



x (G°W(R' + T -, R + T -, -e n )g x ,G^\R' - R x + ^, e n ) 
X9» P G° K f(R2 + | , Ri - j, -e n )g Ka G% R \R 2 - | , R - e n )) 



x / ^'gjf 2 A:,(R, k')A^(R, ki)A*,„(R, k 2 )e'( k "^^^^ . (2.32) 



Here II = — jVr- 2eA w /(R) and we have assumed the slow variation of the vector potential. 
c '"""5 G Q a p = G°5 a (3 is diagonal in the spin space, we can express A^'^p = A*^^ g a p 



with 

A*f <*> (R, k; id) = A*?M (Rj k . w) + A *£(*) (R> k . w) j (2 . 33) 

where 

dk f 



A*™(R,k) =T £ / (f^( k " k ') / dre"*'*/ / rfR'rfr' 



en>0' 

x (G^R' + 1 R + J; -(e n + ,))^y(R' - R - J; e n )) 



x/^ e ik "' r 'A:(R,k"), (2.34) 



(27T) 

and 
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A* 9 (R,k) = -iT £ / -^-F(k-k') y rfre- k '- r | j dR'dr' 



e n >0' 



x (G^(R' + 2 ' R + 2 ; " (e " + CJ ))^G r y (R' - 2 ' R - 2 ; €n)) 

r die" 

x[(R' - R) • II] 2 / |^e ik "- r 'A:(R,k") , (2.35) 

bimilariy, ZX J/a/3 - l\ n g aj3 . 

Using the diagrammatic rule mentioned above, we calculate A*j'^ R \ which in general 
has four terms for the mixed s and d wave superconductors, Aj^ ^ = Y^t=i AjJ, ^ . It is 
easy to show that 

A*^(R,k;u;)=T£ / -^- V ^V s G^(-p; -(e n + W ))G°W(p; e n )A s *(R, u ) 



e„>0' 



^^{[ln ^ + ^(i) - + £)] + ^'(- + ^)} 
2 u ttT rv 2 ; v 2 2 ;J 4vrT v 2 2 

xA;(R,w) , (2.36a) 



e n >0 • 



(2 



= , (2.36b) 
A$ (fl) (R,k;u;) = T £ / ^^V^-^)(pt-p 2 )G°^(-p;-(6 n + a;))G W(p; eri ) 



en>0' 



<ip 

(27) 

xA*(R,cj) 

, (2.36c) 



AgT(R, M = TE/ 7^(*2 " - fy)G 0{R) (-p; -(e n + *))G°«>(p; e n ) 



xA*(R, W )(^-@ 

= M^{[ln + - ^(1 + *)] + ^'(- + ^)} 

xAJ(R )W )g4;), (2.36d) 

where 7 is the Euler constant, uo is the cut-off frequency, p s = 1/iyTt s , p\ = 1/27TTi, -^'(x) 
is the derivative of the digamma function ip(x). 

As for the results of A% N {R) = Et=iK' 9 f R) and A*f {R) = EHiKiTl the second 



term in Eq. fl2~27D can be dropped since it gives a very small higher order correction due to 



the fact that the gap function usually has a temporal variation over a time scale very long 
compared to the range of the Green function, that is, u)T\ <C 1. Therefore, the details to 
evaluate them are the same as the static caseJ^M Here we just give the results 
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A^(R,k;u;) = -^^(^[x^A^R; u>) + l X i, 2 (n 2 - IlJ)A*(R;u/) 



1 — ' - [xi^in^ - ry a s (r ; u) ) 



16 VT' 

+Xo, 3 n 2 A^(R;^)] (2.37) 



A*f (ii) (R,k;a;) = -MM{[ X3>0 - p. X4i0 ]A:(R; W ')A.(R; W w )A:(R;a;* 
+IX2,i - ^x^JA^R^OA^R^'OA^R;^") 
+1x2,1 A^(R; wO A S (R; u") A*(R; a/"] 



V^(g - fc 2 )iV(0) 
4(vrT) 2 



{X2,iA:(R; a/)A d (R;a/')A;(R; u/' 



+[X2,i - ^]AS(R;wOA.(R; t OA:(R;w w ) 
+^ X o, 3 A^(R; u/)A d (R; uj'')A* d (R; u'"]} , (2.38) 
where v p is the Fermi velocity, and Xm,m' is a function defined as 

Xm ' m ' = 5 (2n + 1 + p s )™(2n + 1 + Pl ) m ' ' (2 ' 39) 



C. Anomalous part 

The anomalous part contains integrals of the products of the retarded and advanced 
Green functions, and is therefore sensitive to the details of the spectrum. Following Ref.i, 
we summarize here the diagrammatic rule for the evaluation of this part. In each diagram, 
the solid (electron) lines forming the upper part of the diagram correspond to the retarded 
Green function G°( R >(p; e) = [e — £ p +i/2ri]~ 1 for those lines with arrows to the right and to 
G°^ R '(—p; — e) = [— e — £_ p — z/2ti] _1 for those with arrows to the left. The solid lines in the 
lower part of the diagram correspond to the advanced Green function G 0( - A ^(p; e) = [e — £ p — 
i/2-Ti]" 1 for those lines with arrows to the left and to G 0( - A \— p; — e) = [— e — £_ p + i/2ri]~ 1 
for those with arrows to the right. The triangle and the thin wavy line represent the order 
parameter A and the vertex interaction with the electromagnetic field, respectively. The 
dashed line corresponds to the impurity scattering. If the dashed line encompasses an even 
number of A, a factor l/27rr 1 A^(0) should be assigned. If it encompasses an odd number of 
A, a factor l/27rr 2 iV(0) should be assigned. 

As shown in Fig. [3], the staircase which is the summation of the ladder diagrams, has a 
singular value. We denote it by I(u,k), which satisfies a ladder- type equation 

k) = ^^ {1 + / d P G°^(p; e)G°^(p - k; c - k)} 
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2vrr 1 iV(0) L v ' J J (e - £ + ^)(e - w + v F • k - 



{1 + iV(O) / 1 J / , , /(a;, k)} . (2.40) 



27rriiV(0) L v 7 7 - v F ■ k + i 
Under the condition cjti <C 1 and w^fcri < 1, we obtain 

2^X1^(0) (— iu + Dk z )Ti 

where 7J = v 2 F Ti/2 is the difussion constant for the two dimensional systems. In the real 
time and coordinate space, I~ l is proportional to the operator ■§? — DV 2 . It is important to 



note that the denominator of Eq. Q2.41 ) can be very small if ut\ and Vpkri are both small. 



This fact makes it necessary to sum additionally diagrams containing arbitrary number of 
staircases k), separated by parts including A and A*. Under the assumption t s A s <C 1 
and 7~iAi <C 1, we need only be concerned with those diagrams of the order of A 2 . These 
diagrams together lead to the diffusion equation for the vertex parts T + and T~ as shown 
in Fig. |. 

The kernel Qi corresponds to the diagrams shown in Fig. |3] and is given by 



with 



Qi(R;u) = Q[ a) (R; u) + Q? (R; u) + Q[ c) (R; u) (2.42) 



Q?>(fr«>) = -j ^[G^)(p; e )] 2 G°^(-p;- e )G^)(p; e ) 

xA^(R,p)A:(R,p), (2.43) 

Q?<&«) = -^n(0) J ^^)(p;e)G^)(-p;- e )G^)(p; e )A a/ (R,p) 



x 



/ 7^G !0(ii) (P , ;e)G' ^(-p';-6)G W(p';6)A:(R,p') , (2.44) 



^(h^) = -^/^[^(pV)]^)(pV) 

x / ^[G ^( P ;6)] 2 G°^(- P ;- e )A^(R !P )A:(R, P ) . (2.45) 

Here for simplicity of notation, the vortex renormalization factor is included in the order 
parameter. For the ci-wave component, it is unrenormalized, i.e., A^ = A^; while for the 
s-wave component, A s = r)^ R ' A ^A s or A s depends on whether the vertex connects only 
retarded and only advanced Green functions or it connects a retarded and an advanced 
Green function. 

In view of the specific form of the s-wave and ci-wave pairing states, we see that the 
contribution from the s-wave and ci-wave component is decoupled. Therefore, we obtain 
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Qi(R; u) = Qi,,(R; w) + Q M (R; a;) 



As(R . ^ A:(R . w) _ glf^ A^R; a/)A^(R; w ) . (2.46) 



We can find Q 2 from Q\ by merely replacing all explicit z's by —z's 

Q 2 (R; U ) = ^d^P- A^R; u')A* s (R; u) + S^j^ A rf (R; c/)A^(R; w ) (2.47) 

r s V 2ri / 

The diagram shown in Fig. ^ leads to 

Q 3 (R; w) = (R; o>) + Q? (R; w) + Q 3 C) (R; w) , (2.48) 

with 

QP(BW) = J ^G ^(p; e )G^)(-p;- e )G^)(p; e )G°^(-p;- e ) 

xA^(R,p)A:(R,p), (2.49) 

Q? (Rj") = ^r^y / 7 ^G^ ) (p;e)G^)(-p;- e )G^)(-p;- e )A^(R,p) 

x / ^G ^(p'; e )G ^(p'; e )G ^(-p';- e )A:(R,p') , (2.50) 

SSW) = ^7o) / 7^^)(p; e )G^)(-p;- e )G^(p; e )A.,(R,p) 

x / ^G ^(-p';- e )G ^(p'; e )G^)(-p';- e )A:(R !P ') . (2.51) 

The algebra gives 

Q 3 (R = 2 ^ffi T ' A.(R; Wx)A:(R; a*) + -^^- ^(R; ^)A*(R; a*) . (2.52) 



From the results of Qi, 2 , 3 given by Eqs. fl2.46| ), ( [2.47D , and ( |2.52| ), it is not difficult to prove 
the relation 

Q 3 (R; uj) = -(Qi(R; w) + Q 2 (R; w)) . (2.53) 
In addition, the other three separate terms follows: 

Si(R;u,) = - / ^G ^(p;e)G ^(-p;- e )G ^(p; e ) 

x(tanh^-tanh^^)A Wl (R,p)A: 2 (R,p) , 

= -^cosh- 2 ^= f T ] N f ) A s (R;a;i)^ 2 A;(R;^ 2 ) 
2i 21 e + l/T s 

+ 2(7+ i I/ ( 2°! 1 ) A " (R; ^ V2AS(R; ^ } ' (2 - 54) 
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S 2 (R;c) = -/ ^G ^(p;6)G ^(-p;-6)G^)(p;6) 
e e — ui 



x(tanh— -tanh^^)A wl (R,p)A: 2 (R,p) , 

2T C ° Sh ~ 2 2^ [ T^f As(R; ^ 

+ 9 r iiV / ( 9 0) A ( R ' ^)^A*(R; u 2 ) , (2.55) 



and 



S 3 (R;") = / ^(tanh^-tanhi^)G ^(p; e )G ^(p; e )e^(R) 

= ntT ]W® cosh" 2 ^(-^)e^(R) . (2.56) 

Here we have approximated 

, e e — u; a; , _ 2 e 

tanh — - — tanh — — — — - cosh — - , (2.57) 
2T 2T 2T 2T y ' 

when < T. 

From the results for Q's and S"s, we obtain the diffusion equation for T + 



4Tti 2T l e 2 + r- 2 dt 2[e 2 + (2t 1 )- 2 ] dt dt 

^s 1 , A «9A^ A *<9A S , (2T!)- 1 <9A^ A *dA d . 



e 2 + r s - 2V 6 s <9t 7 2[e 2 + (2r x )- 2 ] v a dt d 9t 

-{^l A »l 2 + 2l Jltn)^ A ' ?}{r+ + T ~ h (2 ' 58a) 

where A Sj< 2 and </? are functions of R and i. and are funtions of R, t, and e. 
Similary, the diffusion equation for F~ is found to be 

^dt " ^ ^ ^ 4T^ C ° Sh 2T<? + ^ ^ + 2 [g 2 + (2r 1 )- 2 ] ft " ^ 

r; 1 (A OA: dA s (2TQ- 1 9AS 9A, 



e 2 + r- 21 s dt s dt ' 2[e 2 + (2n)- 2 ] 1 d dt d dt 

-^' A -' 2 + 2[^X-^ A -' 2 » r+ + r '»- < 2 ' 58b » 
These two diffusion equations can be rewritten as 

d_ _ _ = *cosh~ 2 (^) e d\A s \* d\A d \\ 

y 8t n ' 2TT! e 2 + r- 2 ft 2[e 2 + (2r x )- 2 ] ft 

(2.59a) 
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2[e 2 + (2r,)- 2 ] 1 "at " 8i ' 9i 

-<^T^I A «I 2 + e . + ff)-. !^!'}^ + r ") ■ < 2 - 59b » 

With the results of T ± , the anomalous part represented in Fig. |5| is given by 

e = ^{-/&/^v(k- P )G°^(-p^ 

-/ de J ^V(k-p)G ^(-p;- e )G ^(-p;- e )G ^(p; e )A: i (R,p)rj 2 (R,e)^ 

= -^^a:(R;^) /^i^i- (R, e) - ^r- 2 (R, e)] 

- ^)A*(R;^) / de[ ^^ r+ 2 (R,e) 
e + Z(2ri Sr j2 (R, e) ]. (2.60) 



e+(2r!)- 2 



where is a function of R, u, and k. 



D. TDGL equations for the order parameters 



From Eqs. (|2.36 ), ( 2.37|) , and ( |2.38|) , and fl2.60|) , by performing the inverse Fourier trans- 
form and comparing both sides of the gap function for k-independent terms and the terms 
proportional to k 2 x — A; 2 , we obtain the coupled TDGL equations for the order parameter 
components: 

_X^p) +2 $ s ( R) t)A:(R,t) = 
2a s A:(R,t) + ) 2 [I X2il n 2 A:(R,t) + | Xl , 2 (II 2 - II 2 )A*(R,t)] 



+ 



•■■kT 
Iff) 



{2( X 3,o-p s X4,o)A:(R,t)|A s (R,t)p 



+ 2(x 2 ,i - f x 2 ,2)|A,(R,t)| 2 A:(R,t) + X 2,iA5 2 (R,t)A fl (R,*)} 



(2.61) 



-^^f^ + 2$,(R,t)A S (R,t) 



Id 

a d A* d (R, t) + (^) 2 [ X0i3 n 2 A*(R, t) + xiAK - n 2 ) a:(r, t)) 



+ 



TIT 1 ) 



{(Ixoj A* d (R,t)\A d (R, t)\'< 



+ 2( X 2,i - f X2, 2 )|A s (R,t)| 2 AS(R,t) + x 2 ,iA: 2 (R,t)A,(R,t)} . 
Here ~f s d are two relaxation rates defined by 



(2.62) 



7. 



-i 



2vrT rv 2 2' 



(2.63) 
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and 



= iSr*<5 + 7 > • (2 - 64) 



Two quantities $ s ^ are given by 



iri /" , r e — it 1 , e + ir 1 



/i s (R, t) = --± de[ % r+(R, t, e) - _-^r"(R, t, e)] , (2.65) 

A^OM = -^/ d e [ ^- + i ^p I^(R,t,e) - ^ ( ( ^ ) ) ^ r-(R,t,e)] ■ (2.66) 
Finally, the parameters a s and are given by 

«, = -[ln^ + ^)-^ + f)], (2.67) 



and 



« rf = -[ln^? + ^)-^ + f)]- (2.68) 



T csQ and T c d are the critical temperatures of a clean superconductor, which are determined 
by 

N{0)V s \n(2e^ D /nT cs0 ) = 1 , (2.69) 

and 

[N(0)V d /2] ln(2e^/vrT cd0 ) = 1 , (2.70) 

with 7 the Euler constant and u>£> the cut-off frequency. In the presence of impurity scat- 
terings, two transition temperatures are determined by the conditions a s (T cs ) = and 
&d{T c d) — 0. It is very clear that the transition temperature T cs for s-wave order parameter 
can only be affected by the magnetic impurity scattering while the transition temperature 
for rf-wave order parameter is dominantly affected by the nonmagnetic scattering. The crit- 
ical temperature of the superconductor is defined by T c = max{T cs , T C( i\. We estimate that 
as long as the <i-channel interaction Vd is larger than about three times of the s-channel 
interaction V s , the pure c?-wave state is stable in the bulk systems without perturbations. 
The phase diagram of such a system in the absence of external fields and impurities has 
been previously studied in Ref.EJ. 

By introducing a formal free energy density 

/(R, i) = 2a s | A S (R, t)| 2 + « d | A d (R, i)| 2 + ^ 

nil 4 

+^xi,2[n:A s (R, t)n x A*(R, t) - ii;a s (r, t)n,A*(R, t) + c.c.]} 

+(^) 2 {(X3,o -p,X4,o)|A s (R,t)| 4 + jj( X0i3 - ^i X04 )| Ad (R,t)| 4 

+(2x2,1 -piX2,2)|A d (R,t)| 2 |A s (R,t)| 2 + i X2il [A* 2 (R,t)A 2 (R,t) 
+Af(R,t)A 2 (R,i)]} " (2.71) 
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the TDGL equations (|2.61|) and ( |2.62| ) can be written in a compact way 

iaA ^ + 2*.(R J t)A:(R,t) = ^ I (2.72a) 



7s dt »v > / ,v > , 



Id dt »v w av w 5Ac 



III. TIME-DEPENDENT CURRENT AND CHARGE DENSITY 

A. Current density 

The expression for current in "imaginary" frequency space is given as 

J(x,r) = - — (V x - V x 0(5G(xr;xV o+ ))| x ^ x - — A(x, r) (5G(xr; xr 0+ )> , (3.1) 

where {5G a p) = (G a p — G° a g) = (SG)5 af 3 with G a p defined by Eq. (|2.10|) and the factor 2 
arises from the spin sum. Using the similar technique for the gap function, we can divide 
the current into the normal and anomalous parts, that is, 

J u (R) = J?(R)+Ji(R,t). (3.2) 

The normal part is given by 

,e7\ 



ml e n >0 



£ J »Wr"[e- ! ( R '- R )- n * + ( r '- r '- Vr A QJ1 (R,r)] 



x 



[ e i(R"-R)-n+(r»-r).v rA * (R> r )] G o(fl) (R // _ r // /2; R / _ r / /2 . _ £r 



^2 * 

x V r {G°^(R + r/2, R' + r'/2; e n )G°^(R" - r"/2, R - r'/2; e„)} r ^ 



'mi' 



X] / dRWdR'tir'^e-^^-^-^+^-^-^A^^R^)] 



e n <0' 

x [ei (R"-R).n+(r"-r).V rA * 2(R) r )] G 0(A) (R /, _ ^ R / _ ^3. ^ 

x V r {G°^(R + r/2, R + r'/2; e n )G°^(R" - r"/2, R - r'/2; e„)} r ^o ■ (3.3) 
The computation of this part is the same as the static caseEHEH and we give as a result 

(R) = ^f^4x2,iA:(R; ^)ITA S (R; atf + ^Xo, 3 A* d (R; ^ 2 )ITA,(R; u x ) 
m{Tcl ) 2 4 

+ ^Xi, 2 [A:(R; cu 2 )fKA d (R; u x ) + A*(R; w 2 )n;A s (R; u x )]e x 

-A:(R; ^ 2 )n;A d (R; Wl ) + A*(R; a; 2 )n;A s (R; Wl )]e y ]} + C.C. . (3.4) 

Here e XjV is the unit vector along the x(?/)-direction. 

The anomalous part is represented by the diagram shown in Fig. || The contribution 
from the first term is 
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2C ! f 6 6 — Ul f 6 

R ^ = "^m^Ivri J rfe ( tanh 2T ~ tanh I ^'["^M + ^ A ^( x ") ' P*"l 

xV r [G 0(R) (R + r/2, x"; e)G 0(A) (x", R - r/2; e)] r ^ 



x[-e V „(x") + -A„(x") ■ P2 ] Vr [e->"'< R+ '/ 2 - x "> e ,K '< II "- R+ " 2 >; 



??? 



2 ^ 1 /rf e (tanh^-tanh^)/dx''/^G°^(p;6)G^)(p; e ) 



mi' Ani J y 2T 2T J J (2vr 



Q-T X , 
2ie 

4ie 



xp[-e^(R) + — A W (R) ■ p] 
m 

= -a[-^A w (R)] , (3.5) 

where a = N(0)e 2 VpTi = 2N(0)e 2 D is the normal-state conductivity. Here we have used 
the integral 

de^- cosh" 2 -J- = 2 . (3.6) 

2T 2T v ; 

Similarly, the contribution of the second term is given by 

x V r [G 0(i?) (R + r/2, x"; e)G 0(A) (x", R - r/2; e)] r ^ 
V J deT+(R,e) 

v/de[r+(R,e)+rj(R,e)]. (3.7) 

By performing an inverse Fourier transform to Eqs. Q3.4| ), ( |3.5| ), and fl3.7| ), we obtain the 
current in real time 

J(R,t) = J„(R,t) + J fl (R,t) • (3.8) 
Here the normal state current is given by 

J„(R, t) = -a[Vm, t) + dM ^ t] ] (3.9) 

where 

<p(R, t) = -^J de[T + (R, t, e) + r~(R, t, e)] (3.10) 

can be considered as the effective electro-chemical potential for quasi-particles. The super- 
current is given by 

J.(R,t) = SE J N J® {^X 2 ,iA:(R, t)IT A.(R, t) + ^o, 3 AS(R,t)n*A d (R,t) 
m(7cl j 2 4 

+±Xi,2 [A: (R, t)n:.A d (R, t) + A* d (R, t)U* x A s (R, t)]e x 
-A* S (R, t)n;A d (R, t) + A*(R, t)n* A S (R, t)]e y ]} + C.C. 

iV(o)*/0M) 

" ~4 6A~ • 
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B. Charge density 



The charge density in the "imaginary" time space is defined by 

p(x, r) = -2e(G(xr, xr +0 )) . (3.12) 
After the analytical continuation, we have 

p u (x) = -2eT^G £ie _ a ,(x,x) 

e 

= /^(R)+p A (R), (3-13) 

with 

(R) = ( _2e )^~~ / rfetanh ^ / ^i(-ey w (xi) 

x [G 0(K) (x, Xl ; e + cu)G 0(fl) (x 1; x; e) - G 0(A) (x, x i; e)G 0(A) (x 1; x; e - u) 

= -e 2 iV(0)^(R) / di{[V— + 5(0} + [-P^ + 6(0}} 

= -2JV(0)eV(R) , (3.14) 



and 



p A (R) = -(A)— / rfe(tanh - tanh ^— ^) / dx"(-e^(x") + eA w (x") • p, 
mi Am J 21 21 J 

x [G 0{R) (K + r/2, x"; e)G 0(A) (x", R - r/2; e)] r _> 
= (-2e)i/i(-r:(R, £ )) / * G0W(e )P )G°W( C ,p) 



zA(0)er! 
2 

,2; 



(27 

| C /e[r+(R, e ) + r-(R,e)] 
2e'A(0)^(R) . (3.15) 

After an inverse Fourier transform, we have the charge density in real time space 

p(R,t) = 2e 2 A(0)[<p(R,t) - cp(R,t)) . (3.16) 

From Eqs. (|3.8|) , (|3.16| ), and (|2.59|) follows the continuity equation V • J + dp/dt = 0. 



IV. DISCUSSIONS AND SUMMARY 



Combined with the Maxwell equations, which couple A and <p with J and p, Eqs. Q2.72| ), 
( P-8|) , ( ft-16 ), together with Eqs. Q2.59 ) constitute a complete set of coupled time-dependent 



Ginzburg-Landau equations, which are our main results. Several features of the above results 
deserve special attention: It is well known that depairing of s-wave superconductors are due 
only to magnetic impurities. However, nonmagnetic impurities can have direct depairing 
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effects on unconventional ci-wave pairing state. Similarly, the relaxation of the s-wave order 
parameter is influenced only by magnetic impurities. Therefore, the magnetic impurities 
as pair-breakers are essential in the derivation of the corresponding TDGL equations for 
conventional s-wave superconductors.0 However, nonmagnetic impurities acting as depairing 
centers can directly affect the relaxation of the ci-wave order parameter. Interestingly, for a 
mixed d- and s-wave symmetry superconductor with a high concentration of magnetic and 
nonmagnetic impurities such that r{T c <C 1 and t s T c < 1, we have 7^" « n, 77 1 ~ r s . In 
this limit, the TDGL equations for the order parameters, Eqs. (|2.72a| ) and ( 2.72b ), become 



-r s [| + 2^(R,t)]A:(R,t) = , (4.1) 



- Tl [| + 2*ey5(R, t)) A* d (R, t) = , (4.2) 

where the coefficients a s ,ad, and Xm,n can also be simplified, but are not explicitly given 
here. 

These set of TDGL equations valid under the strong gaplessness conditions are similar in 
form to that postulated phenomenologicallyi except that the relaxation parameters obtained 
here are 7 S (= r^ 1 ) and 7^ (= rf 1 ) and the usual scalar potential <p is replaced by the electro- 
chemical potential (p. Therefore, the phenomenological TDGL equations are at most valid 
when the superconductor is very dirty with also a high concentration of magnetic impurities. 
If the superconductor is doped only with high density of nonmagnetic impurities (r s ^> 7~i), 
the TDGL equation ( j2.72b| ) for ci-wave component is reduced to Eq. ( [4.2| ) while the relaxation 
parameter involved in the equation for s-wave component becomes 77 1 ~ tt/AT c . In this 
case, 7 S <C 7^ and the TDGL equations for both components are quite asymmetric. Of 
particular interest, if T cs < T < T c d, due to a mixed gradient coupling of the s- and d- 
wave components, the s-wave order parameter with four-lobe structure is induced near the 
ci-wave vortex core, and the overall structure of an individual vortex is fourfold symmetric. 
Numerical simulation,! where the same relaxation rate (7^ = 7d = 7) was assumed for 
two components, showed an intrinsic contribution to the Hall angle caused by the lack 
of complete rotational symmetry in ci-wave superconductivity. In the case r s ^> n, we 
could have 7 S <C 7^ and the ci-wave order parameter relaxes much faster than the s-wave 
component. Under this condition we expect that the s-wave component will not be able 
to follow the motion of the ci-wave vortex and novel phenomenon may appear in the flux 
dynamics. Even when t{F c and t s T c are both small, the condition 7^ = 7,, used in Ref.i 
would require the assumption that the non-spin-flip interaction JJ\ = 0, which as judged 
from the studies on conventional s-wave superconductors,^ may well be not justifiable. 

In summary, we have derived the TDGL equations for superconductors with mixed ci- 
wave and s-wave symmetry assuming a weak gapless condition for both types of order 
parameters. From this derivation, the unknown coefficients for the TDGL equations postu- 
lated phenomenologically have been ascertained. This set of TDGL equations can be used as 
the starting point for the study of the vortex dynamics in superconductors with the mixed d- 
and s-wave symmetry, or even extended to study other transport coefficients. In particular, 
the issue of how the dynamic properties of vortices are influenced by the admixture of an 
induced s-wave component with the dominant rf-wave component of the order parameter as 
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well as their different responses to the impurity scatterings can be studied systematically. 
The TDGL equations for <i-wave superconductors with on-site s-wave repulsive interaction 
can be similarly obtained by using the Pade approximation,!! and we find that the main 
conclusion of the present paper still remains unchanged. This result together with a de- 
tailed derivation will be presented elsewhere. Finally, we remark that the present derivation 
has not included the effects of electron-electron (actually hole-hole), electron-phonon, and 
electron- "any magnetic excitation" scatterings, which might be more important in high-T c 
superconductors than in conventional low-T c superconductors. Whereas such inelastic scat- 
terings are far from being easy to incorporate within the present framework, we think that 
their dominant qualitative and perhaps semi-quantitative effects can be taken into account 
phenomenologically by adding a term 1/te to the diffusion operator d/dt — DV 2 , where te 
stands for an inelastic relaxation time (assuming that the weak gaplessness conditions are 
still satisfied). Consistent with such an approach one should regard t% and r s as effective, 
including also some effects of the inelastic scatterings. 
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FIGURES 



FIG. 1. Ladder-type diagram leading to I(u,"k). The momenta and frequencies for the solid 
(electron) lines in the upper part are p and e, and those for the solid lines in the lower part are 
p — k and e — u. 

FIG. 2. Impurity-averaged diagrams leading to the diffusion equation for r + . The thick wavy 
lines correspond to k) shown in Fig. |l[ The thin wavy line corresponds to the vertex interaction 
with the electromagnetic field. The triangle represents the order parameter. 

FIG. 3. Impurity-averaged diagrams for kernel Q\. A and A* are both the vertex-renormalized 
order parameters in the upper part. 

FIG. 4. Impurity-averaged diagrams for kernel Q3. A is the vertex-renormalized order param- 
eter in the upper part and A* the vertex-renormalized order parameter in the lower part. 

FIG. 5. Impurity-averaged diagrams leading to the anomalous part in the TDGL equation for 
order parameter. are given by the type of diagrams shown in Fig. 

FIG. 6. Impurity-averaged diagrams leading to the anomalous current density. 
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